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Abstract 

We consider the problem of enumerating periodic cx-juggling sequences of 
length n for multiplex juggling, where a is the initial state (or landing schedule) 
of the balls. We first show that this problem is equivalent to choosing I's in 
a specified matrix to guarantee certain column and row sums, and then using 
this matrix, derive a recursion. This work is a generalization of earlier work of 
Fan Chung and Ron Graham. 



1 Introduction 

Starting about 20 years ago, there has been increasing activity by discrete mathe- 
maticians and (mathematically inclined) jugglers in developing and exploring ways of 
representing various possible juggling patterns numerically (e.g., see [H El [3l IH [3 El 
[Tot IT2l HI]). Perhaps the most prominent of these is the idea of a juggling sequence 
(or "siteswap", as it is often referred to in the juggling literature). The idea behind 
this approach is the following. For a given sequence T = {ti,t2, . . . ,t„) of nonnega- 
tive integers, we associate a (possible) periodic juggling pattern in which at time i, 
a ball is thrown so that it comes down at time i + ti. This is to be true for each 
i,l < i < n. Because we assume this is to be repeated indefinitely with period n, 
then in general, for each i and each k > 0, a ball thrown at time i + kn will come 
down at time i + ti + kn. The usual assumption made for a sequence T to be a valid 
juggling sequence is that at no time do two balls come down at the same time. This 
assumption results in many consequences, e.g., all of the quantities i + ti (mod n) 
must be distinct, the number of balls in the pattern is the average (1/"^,) Ylk=i ^k, and 
the number of juggling sequences with period n having fewer than b balls is 6" (see 
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An important object for understanding the relationships and transitions between 
various jugghng sequences is the concept of a state diagram, developed independently 
(and almost simultaneously) by Jack Boyce and Allen Knutson [9J . This is a directed 
graph where each vertex is called a state or landing schedule, a 0-1 vector indicating 
when the balls that are currently in the air will land, and edges represent possible 
transitions between states. The vertex and edge sets for the state diagram can be 
defined as follows: 

V = {{ai,a2,as, ...): UiE {0,1}, Y,iO'i = b}, 

E = {{ai,a2,as, ...)-^{bi,b2,h,...) : ai<bi_i foi i = 2, 3,...}. 

More specifically, each juggling sequence T is associated with a state a = = 
(fXi, (J2, . . . , CT/i, . . .) which can be found by imagining that the sequence has been 
executed infinitely often in the past, with a final throw t„ being made at time 0. 
Then is 1 if and only if there is some ball still in the air at time that will land at 
time i. In this case we say that T is a a-juggling sequence. 

If we are now going to throw one more ball at time 1, transitioning to a (possibly) 
new juggling state a' , then we are restricted to throwing it so that it lands at some 
time i which has aj = 0. The new state a' = {a[,a2, ■ ■ ■) then has a[ = cr^+i for 
k>l,ky^j — 1 and o'j_i = 1. The preceding remarks assume that ai = 1. If o"i = 0, 
so that there is no ball available to be thrown at time 1, then a "no-throw" occurs, 
and the new state vector a' satisfies = ak+i for all > 1. These give the two basic 
transitions that can occur in the state diagram. With this interpretation, it is easy 
to see that a juggling sequence of period n exactly corresponds to a walk of length n 
in the state diagram. 

In [5], the problem of enumerating cr-juggling sequences of period n was studied, 
which by the above comments is equivalent to counting the number of directed closed 
walks of length n starting at a in the state diagram. In the same paper, the related 
problem of counting the number of "primitive" closed walks of length n was also 
solved. These are walks in which the starting state a is visited only at the beginning 
and the end of the walk. 

A particular unsolved problem mentioned in [5] was that of extending the analysis 
to the much more complex situation of multiplex juggling sequences. In a multiplex 
juggling sequence, for a given parameter m, at each time instance up to m balls can 
be thrown and caught at the same time, where the balls thrown at each time can 
have different landing times. Thus, ordinary juggling sequences correspond to the 
case m = 1. 

As before, we can describe a (multiplex) juggling sequence as a walk in a state 
diagram. Here a state a = (oi, 02, 03, . . .) can again be described as a landing schedule 
where a, are the number of balls currently scheduled to land at time i. We also have 
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a state diagram which has as its vertices all possible states and for edges all ways to 
go from one state to another state (see [ID]). The state diagram is thus a directed 
graph with two important parameters: b, the number of balls that are being juggled, 
and m, the maximum number of balls that can be caught /thrown at any one time. 
The vertex set and edge set are defined as follows: 

V = {{ai,a2,a3, ...): ai e {0,1,. .. ,m},Y,iai = b}, 

E = {(01,02,03, •• •)-^(^i,&2,&3, •••): flj < for i = 2,3, .. .}. 

Since each state will only have finitely many nonzero terms, we will truncate the 
terminal zero portions of the state vectors when convenient. The height of a state a 
will be the largest index i for which Oj > 0, and will be denoted by h{a). A small 
portion of the state diagram when 5 = 3 and m = 2 is shown in Figure [T] 

(l,0,2).i (2,1) O ^ 

(1,2) ;h (1,1,1)0 (0,2,1) 

(0,1,2)^1 (2,0,1) 

Figure 1: A portion of the state diagram when 6 = 3 and m = 2. 



1.1 A bucket approach 

To better follow the analysis (and make practical demonstrations easier) we can rein- 
terpret multiplex juggling by a series of buckets and balls. The buckets will represent 
future landing times for i = 1, 2, 3, . . . and the balls are distributed among these 
buckets. A state vector is then a listing of how many balls are currently in each 
bucket, and m is now the maximum number of balls that can fit inside of a bucket. 
Transitions from state to state happen by having the buckets shift down by one and 
redistributing any balls that were in the bottom bucket. This process is shown in 
Figure |2| 

1.2 Multiplex siteswap notation 

To describe a walk in the state diagram it suffices to know what state we start in 
and how we transition from state to state. In transitioning from state to state the 
important piece of information is what happened to the ball(s) in the bottom bucket. 
This can be described by a multi-set which lists the new location(s) of the ball(s). 
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Initial state 
= (2,1,0,3,2) 



QQ/ 



Buckets advance one 
step freeing balls in 
the bottom bucket 



QQ/ 



QQ/ 



Freed balls are redistributed 
to available buckets to get the 
next state = (2, 0, 3, 2, 0, 1) 



Figure 2: A bucket approach to multiplex juggling. 



We can thus describe our walk by a series of multi-sets (Ti,T2, . . . ,T„) such that 
each set has m elements (when we have fewer than m balls to redistribute we will 
indicate no-throws by 0). These sets are analogous to siteswap notation for juggling. 
In particular, it can be shown that 



[j{^ + T, (mod n)} = {l,...,l,2,...,2,... 



1 

} , and — X = b. 



i=l 



n 



rn times m times 



1=1 xeT, 



In the next section we will combine the idea of this multiplex siteswap notation 
with the buckets. 



2 A matrix interpretation 

One way to use the buckets to find a sequence of length n that starts in state a = 
(oi, a2, . . .) and ends in state j3 = 62, • • •) (if one exists) is to start with the balls 
arranged in the buckets as dictated by a. We then modify the capacities of the 
buckets so that they are (starting at the first bucket) 

m,m, . . . ,m, 61,62, 

V ' 

n times 

Finally, take n steps (such as shown in Figure [2]) being careful not to exceed the 
capacity of any bucket and at the end, we will be forced into state (3. On the other 
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hand, every possible way to start in a and end in /3 in n steps can be done in this 
modified buckets approach. 

Finding all of the walks of length n in the state diagram between a and /3 is 
thus equivalent to finding all of the walks that can be run using this modified bucket 
procedure. This is what we will actually enumerate. We start with the following 
matrix, where h = max{h{a) — n, h{(3)} (this is a 0-1 matrix, where any unspecified 
entries are O's). 



A 



o 
+ 



m columns 



1 1 
1 1 
1 1 



\l 1 



1 

1 1 1 

1 1 1 

1 1 1 



1 1 1 



colu 



m columns 



1 1 
1 1 



111 



1 1 



n blocks each with m columns 



Each block of m columns will correspond to one transition in the state dia- 
gram/buckets procedure. In a block, each column corresponds to a single element 
in the multi-set describing our transition. The first 1 in a column corresponds to a 
no-throw or a throw of height 0, the second 1 corresponds to a throw of height 1, the 
third 1 to a throw of height 2 and so on. So we choose exactly one 1 in each column 
and by reading the blocks of columns we will get the transitions Tj between states. 

Each row will correspond to a bucket, and to incorporate the modified buckets 
approach we specify a row sum for each row. The row sum will be the "unused 
capacity" of the buckets. Beginning at the first row and going down the row sums 
will be 

Oi, m — . . . , m — a„, bi — a„+i, . . . ,bh — an+h- 



m 



As an example, suppose that we wanted a walk of length 4 that starts and ends 
in the state (1,2) in the state diagram shown in Figure [T} Then this corresponds to 
choosing one 1 out of each column in the matrix below on the left so that the row 
sums are as dictated on the side of the matrix. 
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One possible solution is shown in the matrix on the right which corresponds to the 
walk 

(1,2) ^ (2,1) ^ (2,0,1) ^ (0,1,2) ^ (1,2). 

{0,2} {1,3} {3,3} {0,0} 

Here we have indicated the sets Tj under each transition. 

To see why this works we need to show that each multi-set given by a block of 
columns is a valid transition in our modified buckets procedure. Starting with the 
first block, the first row sum should be m — «!, this indicates that we currently have 
an excess capacity of m — ai in the bottom bucket and so we need to make m — «! 
no-throws, i.e., throws of height 0, and by row sum restrictions we are forced to select 
exactly m — ai of the I's on the first row. This accounts for m — ai of the columns 
in the first block. The remaining ai columns must have I's (i.e., balls) distributed 
among the rows (i.e., buckets) which still have extra capacity. Thus the first block 
must give a valid transition in the procedure. After selecting the I's for the first block 
we then update the row sums according to our choices, remove the first block and 
then repeat the same process n — 1 times. 

Conversely, it is easy to check that given the transitions Ti, T2, . . . ,Tn joining a 
to f3, we can find a selection of I's in the matrix A which corresponds to this walk 
and satisfies the row/column sum restrictions. 

Finally, since multi-sets are unordered, our choice of I's in columns is unique up 
to permutation of the columns in a block, so we may assume that the heights of the 
I's in a block are weakly decreasing (as shown in the example above). We now have 
the following general result. 

Lemma 1. Suppose we are given a state diagram with capacity m and states a — 
(01,02,0.3,...) and (3 = (61, • • •); (ind h = max{/i(Q;) — n, /?.(/?)}. Then the 

number of walks of length n starting at a and ending at (3 is equal to the number of 
ways of choosing mn ones out of the matrix A given above such that: 

• Each column sum is 1. 

• The row sums are (in order from first to last) 

m - oi, m - 02, . . . , m - a„, &i - a^+i, . . .,bh- On+h- 

• In each block of m columns the height of selected 1 's is weakly decreasing. 

Remark 1. If for some i, hi < a^+i, then one of the row sums will be negative which 
is impossible and thus we have no solutions for the selection of I's. At the same time 
it is easy to see that there can be no walks in the state diagram of length n joining 
a and (3 by comparing their landing schedules. 
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Remark 2. When m = 1, if we ignore rows with row sum 0, then all the rows and 
column sums will be 1. In this case we can count the number of walks joining two 
states by calculating the permanent of a matrix. This is similar to the approach taken 
by Chung and Graham |5j. 



3 Filling the matrix 



We now count the number of ways to fill the matrix A according to the restrictions 
in Lemma [Tj We will demonstrate the procedure by working through an example, 
namely counting the number of periodic multiplex juggling sequences of length n that 
start and end in state (3) when m = 3. 

The first thing to observe is that when the height of a is small compared to n 
then the row sums have the following form: 



initial 



noise 



terminal 



noise 



string 
of m's 

We can form a recurrence based on this where we fill in the last block and then reduce 
it to the case when there is one fewer m in the middle. 

Without loss of generality we can assume that our noise at the end is a partition 
of the b balls so that no part is larger than m, i.e., we can ignore any rows with row 
sum and by a simple correspondence we can assume that the row sums are weakly 
decreasing. 

For each partition 7 of 6 with each part at most m, let x^{k) be the number of 
ways to fill the matrix A where the row sums are given by 



m 



ai, . . 



1, and so we will have 



m - ah(a) , m, . , m , 7- 

k times 

In our example, there are 3 such partitions, 3 = 2 + 1 = 1 + 1 
the three functions x^i^k), X2,i{k) and xi^i^i{k). 

We now form a system of linear recurrences, x.y{k) = Y2ga^^sxs{k — 1), by exam- 
ining the different ways to fill in the last block of columns. Note that the row sums 
corresponding to the last block are m, 7. The m comes from the no-throws and after 
we have filled in the last block it will become incorporated into the new terminal 
distribution. Thus at each stage we will decrease the number of middle m's by 1. 

For example, if we are considering X2,i then we have the following 6 different ways 
to fill in the last block satisfying the conditions of Lemma [TJ 

01 



3 


1 


1 


1 


3 


111 


3 


1 


1 1 


3 


1 


1 


1 


2 


1 


1 


1 


2 


1 10 


2 


1 


1 1 
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By looking at the new terminal distributions it follows that 

X2,iik) = xsik - 1) + 4x2,1 (/c - 1) + - 1). 

Similar analysis shows that 

x^lk) = 2x3(fc — 1) + 2x2,i(fc — 1) and 
xi,i,iik) = X3(A; - 1) + 3x2,i(A; - 1) + 4a;i,i,i(A; - 1). 

We can rewrite this in matrix form as 

xsik) \ / 2 2 \ / Xsik-1) 
X2,i{k) =14 1 a;2,i(A;-l) 
xi,i,i(A;) / V 1 3 4 / \ xi,i,i(A;-l) 

In general, we will have that [x-y{k)) = A(^x-y{k — 1)) where A^^s = cl^,s as given 
above. Note that the matrix A will be independent of our choice of a and (3, and 
depends only on b and m. Varying a will change the initial conditions and varying j3 
will change which of the x^ we are interested in for our recurrence. 

For our problem, we are interested in x^ as our terminal state is (3) (i.e., corre- 
sponding to the partition 3). We want to transform our first-order system of linear 
recurrences into a single recurrence for x^. Manipulating the recurrences in ([T]) it can 
be shown that 

X3{k + 3) = 10x3(A; + 2) - 27x3(A; + 1) + 20x3(A;). (2) 

This gives a recurrence for the number of periodic juggling sequence of length k for 
k sufficiently large (i.e., we have shifted our sequences x^{k) past the initial noise; 
when counting periodic juggling sequences we need to remember to account for this 
shift). 

Remark 3. The characteristic polynomial of the matrix in ([T]) is x^ — 10x^-1- 27a; — 20, 
the same coefficients as in (0. This is a consequence of the Cayley-Hamilton Theorem. 
Namely if the characteristic polynomial of the matrix A is x^ + qix^~^ + ■ ■ ■ + qr then 

(0) = {A'' + q,A'-^ + --- + qrl){x^ik)) 

= A' (x^ik)) + qiA''"^ (x^ik)) + ■ ■ ■ + qr (x^ik)) 

= {x^{k + r)) +qi{x^{k + {r - 1))) + ■ ■ ■ + g^(x^(/c)) 

= (x^{k + r) + qiXy{k + {r - 1)) -\ + qrX^{k)) . 

The recursion that the characteristic polynomial of A gives is universal for a fixed b 
and m in the following sense. Let a{n) be the number of walks of length n joining 
state a to state (3. Then for n sufficiently large, a{n) satisfies the recursion given by 
the characteristic polynomial of A (which is independent of a and jJ). 
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It remains to calculate enough initial terms to begin the recursion. Using Lemma[T] 
this can be handled by brute force or some careful case analysis. In general, we will 
need to calculate the first h{a) + r — 1 terms where r is the number of partitions of 
b with each part at most m. The first h{a) — 1 terms are to handle the initial noise 
caused by a and the next r terms are to help start the recursion. A calculation shows 
that the sequence counting the number of periodic juggling sequences that start and 
end at (3) starts 1,4,20. Now applying the recursion we get the sequence 

1, 4, 20, 112, 660, 3976, 24180, 147648, 903140, .... 



With the recursion and the initial values it is then a simple matter to derive a gen- 
erating function for the number of juggling sequences of period n. The generating 
function for this series and several others are given in Table 1. 



State 


m 


Initial terms 


Generating Function 


(2) 


2 


1, 3, 10, 35, 125, 450, 1625, 5875, . . . 


X — 2x^ 
1 — 5a; + 5x'^ 


(1,1) 


2 


1, 3, 11, 40, 145, 525, 1900, 6875, . . . 


') 

X — 2x + X 
1 — 5x + 5^2 


(2,1) 


2 


1,4, 22, 124, 706, 4036, 23110, 132412, . . . 


X - 4:x'^ + 3x^ 
1 - 8x + 13x^ 


(1,1,1) 


2 


1, 3, 18, 105, 606, 3483, 19986, 114609, . . . 


X — 5x'^ + 7x'^ 
1 - 8x + 13x^ 


(2,2) 


2 


1, 3, 21, 162, 1305, 10719, 88830, 739179, . . . 


x-llx' + 33x^ - 27x'' 
1 - 14x + 54x^ - 57x^ 


(3) 


3 


1,4, 20, 112, 660, 3976, 24180, 147648, . . . 


X — Qx^ + Ix'^ 
1 - lOx + 27x^ - 20x^ 


(2,1) 


3 


1,5, 30, 182, 1110, 6786, 41530, 254278, . . . 


X — hx^ + 7x'^ — 3x'^ 
1 - lOx + 21x^ - 20x^ 



Table 1: Number of periodic sequences of length n for some starting states. 



Remark 4. Instead of calculating the first h{a)+r — 1 terms we can instead calculate 
the first h{a) — 1 terms and then compute the r values for the x^(0). In at least one 
case this is easier, namely when the state is (6). Then it is easy to check that x^{fS) = 1 
for all partitions 7 and quickly bootstrap our sequence. 

3.1 Calculating recursion coefficients 

In this section we show how to quickly compute the recursion coefficients 5, which 
can then (by taking determinants) quickly give us the recursion relationship for the 
juggling sequences. 
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Lemma 2. Let 7 and S be partitions of b with no part more than m, and suppose 
that 7 = (71, 72, ... , 7m) and 6 = {61,62, ■■■ , 6m) where 7i is the number of parts of 
the partition 7 of size i, similarly for 6i. Then 



a. 



'7, 




(7i H h 7m) + 1 - {6i+i H V6m) 

6, 



(3) 



Proof. To see ([3]) recall that when filling in the last set of columns of the matrix there 
is the partition 7 and an additional row with row sum m. We then select m ones such 
that (1) the height of the ones are weakly decreasing and (2) no row sum is violated. 

This process can be made equivalent to taking the partition 7 and adding one 
part of size m (to form a new partition 7' = (71, . . . , 7m-i, 7m + 1)) then reducing by 
a total of m some part(s) of 7'. 

The coefficient a^^5 is then the total number of ways that reducing some part(s) 
of 7' by a total of m will result in the partition 6. We can however work backwards, 
namely if we want a desired partition from our reduction we can "insert" the desired 
partition into 7' (i.e., associate each part of 6 with some part of 7' which is at least as 
large) and then finding the difference between the insertion and 7' (which difference 
will sum to m) gives the reduction to use. 

Thus a^^s is the total number of ways that we can insert 6 into 7'. We now 
enumerate by inserting the partition 6 backwards, namely we insert the largest parts 
first and then work down. First note that here are 6m parts of size m to insert and 
they can be positioned into any of the 7m + 1 parts of size m, which can be done in 



ways. There are then 5m- 1 parts of size m — 1 to insert and they can be positioned 
in any of the 7m- 1 + 7m + 1 — parts of size at least m — 1 that have not yet been 
used, which can be done in 



ways. This process then continues, so that in general there will be 6i parts of size i to 
insert and they can be positioned in any of the (7i + ■ ■ ■ + 7m) + 1 — (5i+i + ■ ■ ■ + 6m) 
parts of size at least i that have not yet been used, which can be done in 





(7i H h 7m) + 1 - (5i+i H 6m) 

6^ 



ways. Putting it all together then gives the result. 



□ 
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As a check for, suppose m = 2, 7 = (a, b) and S = (a — 2c, 6 + c), then ([s]) gives 



a + b) + l-{b + c) 
a — 2c 



6+1 
b + c 



a + 1 — c 
a — 2c 



b + 1 
b + c 



To get a nonzero coefficient we must have 6+1 > b+c so that c < 1 and a+l—c > a— 2c 
so that c > — 1. Putting these in and simphfying we are left with 



|a(a — 1) 
(a + l)(6+l^ 
l{b+l)b 




if c 



0; 



if c 
if c= -1; 
otherwise. 



This can easily be verified by hand. 



4 Remarks 

4.1 Primitive juggling sequences 

We have demonstrated a way to count the number of periodic juggling sequences of 
length n which start and end in a given state a. A related problem is counting the 
number of primitive periodic juggling sequences. These are special periodic juggling 
sequences with the extra condition that the walk does not return to a until the nth 
step. This can be done by the following observation of Chung and Graham [5] (see 
also [3 E]): If a{n) is the number of periodic juggling sequences of length n that start 
and end at a and F{x) = J2n>i O'i^)^"', then b{n) counts the number of primitive 
periodic juggling sequences of length n that start and end at a where 

\ b(n)x = 

n>l ^ ' 

Applying this to the data in Table 1 we get the information about primitive 
juggling sequences given in Table 2. Only one of the sequences shown in Tables 1 or 
2 was previously listed in [TT]. 

A related open problem is to find the number of prime juggling sequences which 
start and end in a given state a. A prime juggling sequence corresponds to a simple 
cycle in the state diagram, i.e., it never visits any vertex more than once. Note that 
for a primitive juggling sequence we are allowed to visit vertices other than a as often 
as we want. 
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State 


m 


Initial terms 


Generating Function 


(2) 


2 


1,2,5,14,41,122,365,1094,... 


X — 2x^ 

"1 A ^ 1 o ^2 

i — 4a; + ox 


(1,1) 


2 


1,2,6,17, 48,135,379,1063,... 


1 - 4x + 3a;^ + x^ 


(2,1) 


2 


1,3,15,75,381,1947, 9975,... 


a; — 4x + ox 
1 - 7x + 9x^ + 3x^ 


(1,1,1) 


2 


1,2,13,68,358,1871,9757,... 


X — ox -\- 1 X 
1 - 7x + 8x^ + 7x^ 


(2,2) 


2 


1,2,16,119,934,7463,60145,... 


X - llx^ + 33a;^ - 27a;* 


1 - 13a; + 43x^ - 24a;^ - 27a;* 


(3) 


3 


1,3,13,67,369,2083,11869,... 


X — 6x^ + 7x'^ 
1 - 9a; + 21x^ - 13x^ 


(2,1) 


3 


1,4,21,111,592,3171,17021,... 


X — 5x^ + 7x'^ — 3x* 
1 - 9x + 22x2 - 13x^ - 3x* 



Table 2: Number of primitive sequences of length n for some starting states. 



4.2 Further directions 

One implicit assumption that we have made is that the height the balls can be thrown 
to in the juggling sequence is essentially limited only by the period. This might be 
unrealistic when trying to implement the procedure for an actual juggler. In this case 
we would like to add an additional parameter which is the maximum height a ball can 
be thrown. While it is not difficult to adopt the matrix A to handle this additional 
constraint, our recursion method will no longer work. For this setting, the simplest 
method might be to find the adjacency matrix of the (now finite) state diagram and 
take powers to calculate the number of walks. 

Another implicit assumption that we have made is that the balls are identical, 
it is easy to imagine that the balls are distinct and then we can ask given an initial 
placement of balls and a final placement of balls how many walks in the state diagram 
are there. This problem is beyond the scope of the methods given here. However, 
Stadler [13] has had some success in this direction (using different methods than the 
ones presented here), he was able to derive an expression involving Kostka numbers 
enumerating the number of such sequences, as well as several other related sequences. 

It would be interesting to know for each m, n and b, which states a have the 
largest number of (primitive) cr-juggling sequences of length n. When m = 1, then 
it would seem that the so-called ground state (1,1,..., 1) does. However, for larger 
values of m, it is not so clear what to guess. 

As can be seen there are still many interesting open problems concerning the 
enumeration of multiplex juggling sequences. 
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